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Abstract 

We address an optimal transmit power allocation problem that minimizes the outage probability of a secondary 
user (SU) who is allowed to coexist with a primary user (PU) in a narrowband spectrum sharing cognitive radio 
network, under a long term average transmit power constraint at the secondary transmitter (SU-TX) and an average 
interference power constraint at the primary receiver (PU-RX), with quantized channel state information (CSI) 
^ ' (including both the channels from SU-TX to SU-RX, denoted as gi and the channel from SU-TX to PU-RX, denoted 

as go) at the SU-TX. The optimal quantization regions in the vector channel space is shown to have a 'stepwise' 
structure. With this structure, the above outage minimization problem can be explicitly formulated and solved by 
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, employing the Karush-Kuhn-Tucker (KKT) necessary optimality conditions to obtain a locally optimal quantized 
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power codebook. A low-complexity near-optimal quantized power allocation algorithm is derived for the case of large 



o 

(N 

> 
O 



(N 



number of feedback bits. More interestingly, we show that as the number of partition regions approaches infinity, 
the length of interval between any two adjacent quantization thresholds on the go axis is asymptotically equal 
when the average interference power constraint is active. Similarly, we show that when the average interference 
' power constraint is inactive, the ratio between any two adjacent quantization thresholds on the gi axis becomes 

rS , 

H asymptotically identical. Using these results, an explicit expression for the asymptotic SU outage probability at high 

rate quantization (as the number of feedback bits goes to infinity) is also provided, and is shown to approximate the 
optimal outage behavior extremely well for large number of bits of feedback via numerical simulations. Numerical 
results also illustrate that with 6 bits of feedback, the derived algorithms provide SU outage performance very close 
to that with full CSI at the SU-TX. 

I. INTRODUCTION 

Scarcity of available vacant spectrum is limiting the growth of wireless products and services [[T]|. 
Traditional spectrum licensing policy forbids unlicensed users to transmit in order to avoid unfavorable 
interference at the cost of spectral utilization efficiency. This led to the idea of cognitive radio (CR) 
technology, originally introduced by J. Mitola |l2l, which holds tremendous promise to dramatically improve 
the efficiency of spectral utilization. 



The key idea behind CR is that an unlicensed/secondary user (SU) is allowed to communicate over 
a frequency band originally licensed to a primary user (PU), as long as the transmission of SU does 
not generate unfavorable impact on the operation of PU in that band. Effectively, three categories of CR 
network paradigms have been proposed: interweave, overlay, and underlay [3J. In the underlay systems, 
also known as spectrum sharing model, which is the focus of this paper, the SU can transmit even when 
the PU is present, but the transmitted power of SU should be controlled properly so as to ensure that the 
resulting interference does not degrade the received signal quality of PU to an undesirable level H by 
imposing the so called interference temperature Q constraints at PU (average or peak interference power 
(AIP/PIP) constraint) and as well as to enhance the performance of SU transmitter (SU-TX) to SU receiver 
(SU-RX) link. 

Various notions of capacity for wireless channels include ergodic capacity (for delay-insensitive services), 
delay-limited capacity and outage probability (for real-time applications). These information theoretic 
capacity notions constitute important performance measures in analyzing the performance limits of CR 
systems. In [[51, the authors investigated the ergodic capacity of such a dynamic narrowband spectrum 
sharing model under either AIP or PIP constraint at PU receiver (PU-RX) in various fading environments. 
The authors of [6[ extended the work in [[51 to asymmetric fading environments. In [[3, the authors studied 
optimum power allocation for three different capacity notions under both AIP and PIP constraints. In [[U, 
the authors also considered the transmit power constraint at the SU-TX and investigated the optimal power 
allocation strategies to achieve the ergodic capacity and outage capacity of SU under various combinations 
of secondary transmit (peak/average) power constraints and interference (peak/average) constraints. 

Achieving the optimal system performance requires the SU-TX to acquire full channel state information 
(CSI) including the channel information from SU-TX to PU-RX and the channel information from SU- 
TX to SU-RX. Most of the above results assume perfect knowledge or full CSI, which is very difficult to 
implement in practice, especially the channel information from SU-TX to PU-RX without PU's cooperation. 
A few recent papers have emerged that address this concern by investigating performance analysis with 
various forms of partial CSI at SU-TX, such as noisy CSI and quantized CSI. With assumption of 
perfect knowledge of the CSI from SU-TX to SU-RX channel, [[HI studied the effect of imperfect channel 
information of the SU-TX to PU-RX channel under AIP or PIP constraint by considering the channel 
information from SU-TX to the PU-RX as a noisy estimate of the true CSI and employing the so-called 
'tifr' transmission policy. Another recent work [(91 also considered imperfect CSI of the SU-TX to PU-RX 
channel in the form of noisy channel estimate (a range from near-perfect to seriously flawed estimates) 

2 



and studied the effect of using a midrise uniformly quantized CSI of the SU-TX to PU-RX channel, 
while also assumed the SU-TX had full knowledge of the CSI from SU-TX to SU-RX channel. Recently, 
ifTOl has proposed a practical design paradigm for cognitive beamforming based on finite-rate cooperative 
feedback from the PU-RX to the SU-TX and cooperative feedforward from the SU-TX to the PU-RX. 
A robust cognitive beamforming scheme was also analyzed in [11 J, where full channel information on 
SU-TX to SU-RX channel was assumed, and the imperfect channel information on the SU-TX to PU-RX 
channel was modelled using an uncertainty set. Finally, 111 21 studied the issue of channel quantization for 
resource allocation via the framework of utility maximization in OFDMA based CR networks, but did 
not investigate the joint channel partitioning and rate/power codebook design problem. The absence of a 
rigorous and systematic design methodology for quantized resource allocation algorithms in the context 
of cognitive radio networks motivated our earlier work lOJ, where we addressed an SU ergodic capacity 
maximization problem in a wideband spectrum sharing scenario with quantized information about the 
vector channel space involving the SU-TX to SU-RX channel and the SU-TX to PU-RX channel over all 
bands, under an average transmit power constraint at the SU-TX and an average interference constraint at 
the PU-RX. A slightly different approach was taken in lfT4l . ifTSl where the SU overheard the PU feedback 
link information and used this to obtain information about whether or not the PU is in outage and how 
the SU-TX should control its power to minimize interference on the PU-RX. 

In this paper, we address the problem of minimizing the SU outage probability under an average transmit 
power (ATP) constraint at the SU-TX and an average interference power (AIP) constraint at the PU-RX. 
Similar to |fT3||, we consider an infrastructure-based narrowband spectrum sharing scenario where a SU 
communicates to its base station (SU-BS) on a narrowband channel shared with a PU communicating to 
its receiver PU-RX contained within the primary base station (PU-BS). The key problem is the jointly 
designing the optimal partition regions of the vector channel space (consisting of the SU-TX to SU-RX 
channel (denoted by power gain gi) and the interfering channel between the SU-TX and PU-RX (denoted 
by power gain qq)) and the corresponding optimal power codebook, and is solved offline at a central 
controller called the CR network manager as in iflBll . based on the channel statistics. The CR network 
manager is assumed to be able to obtain the full CSI information of the vector channel space {gi.go) in 
real-time from the SU-BS and PU-BS, respectively, possibly via wired links (similar to backhaul links in 
multicell MIMO networks connecting multiple base stations). This real-time channel realization is then 
assigned to the optimal channel partition and the corresponding partition index is sent to the SU-TX (and 
to the SU-RX for decoding purposes) via a finite-rate feedback link. The SU-TX then uses the power 
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codebook element associated with this index for data transmission. It was shown in [[T3l that without the 
presence of the CR network manager, and thus without the ability to jointly quantize the combined channel 
space, the SU capacity performance is significantly degraded if one carries out separate quantization of gi 
and go. Even if such a CR network manager cannot be implemented in practical cognitive radio networks 
due to resource constraints, the results derived in this paper will serve as a valuable benchmark. Under 
these networking assumptions, we prove a 'stepwise' structure of the optimal channel partition regions, 
which helps us explicitly formulate the outage minimization problem and solve it using the corresponding 
Karush-Kuhn-Tucker (KKT) necessary optimality conditions. As the number of feedback bits go to infinity, 
we show that the power level for the last region approaches zero, allowing us to derive a useful low- 
complexity suboptimal quantized power allocation algorithm called 'ZPiORA' for high rate quantization. 
We also derive some other useful properties related to the channel quantizer structure as the number of 
feedback bits approaches infinity: (a) under an active AIP constraint, the length of interval between any 
two adjacent quantization thresholds on go axis is asymptotically the same, and (b) while when the AIP 
is inactive, the ratio between any two adjacent quantization thresholds on gi axis asymptotically becomes 
identical. Finally, with these properties, we derive explicit expressions for asymptotic (as the number of 
feedback bits increase) behavior of the SU outage probability with quantized power allocation for large 
resolution quantization. Numerical studies illustrate that with only 6 bits of feedback, the designed optimal 
algorithms provide secondary outage probability very close to that achieved by full CSI. With 2-4 bits 
of feedback, ZPiORA provides a comparable performance, thus making it an attractive choice for large 
number of feedback bits case. Numerical studies also show that ZPiORA performs better than two other 
suboptimal algorithms constructed using existing approximations in the literature. Finally, it is also shown 
that the derived asymptotic outage behavior approximates the optimal outage extremely well as the number 
of feedback bits becomes large. 

This paper is organized as follows. Section |Il] introduces the system model and the problem formulation 
based on the full CSI assumption. Section UlI] presents the joint design of the optimal channel partition re- 
gions and an optimal power codebook algorithm. A low-complexity suboptimal quantized power allocation 
strategy is also derived using novel interesting properties of the quantizer structure and optimal quantized 
power codebooks. In Section |IVl the asymptotic behavior of SU outage probability for high resolution 
quantization is investigated. Simulation results are given in Section |Vl followed by concluding remarks in 
Section |VIl 
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II. System Model and Problem Formulation 

We consider an infrastructure-based spectrum sharing network where a SU communication uplink to the 
SU-BS coexists with a PU link (to the PU-BS) within a narrowband channel. Regardless of the on/off status 
of PU, the SU is allowed to access the band which is originally allocated to PU, so long as the impact 
of the transmission of SU does not reduce the received signal quality of PU below a prescribed level. All 
channels here are assumed to be Rayleigh block fading channels. Let gi = and go = |/ioP, denote the 
nonnegative real-valued instantaneous channel power gains for the links from SU-TX to SU-RX and SU- 
TX to PU-RX respectively (where hi and Hq are corresponding complex zero-mean circularly symmetric 
channel amplitude gains). The exponentially distributed channel power gain gi and go, are statistically 
mutually independent and, without loss of generality (w.l.o.g), are assumed to have unity mean. The 
additive noises for each channel are independent Gaussian random variables with, w.l.o.g, zero mean and 
unit variance. For analytical simplicity, the interference from the primary transmitter (PU-TX) to SU-RX 
is neglected following previous work such as flU, ^(in the case where the interference caused by the 
PU-TX at the SU-RX is significant, the SU outage probability results derived in this paper can be taken as 
lower bounds on the actual outage under primary-induced interference). This assumption is justified when 
either the SU is outside PU's transmission range or the SU-RX is equipped with interference cancellation 
capability particularly when the PU signal is strong. 

Given a channel realization {go, gi), let the instantaneous transmit power (with full CSI) at the SU-TX be 
represented by p{go,gi), then the maximum mutual information of the SU for this narrowband spectrum 
sharing system can be expressed as R{gi,p{gQ, gi)) = | log(l + gip{gQ, gi)), where log represents the 
natural logarithm. The outage probability of SU-TX with a pre-specified transmission rate tq, is given as. 
Pout = Pr{R{gi,p{gQ, gi)) < tq}, where Pr{A} indicates the probability of event A occurring. Using 
the interference temperature concept in Q, a common way to protect PU's received signal quality is by 
imposing either an average or a peak interference power (AIP/PIP) constraint at the PU-RX. In [fT6l . it was 
demonstrated that the AIP constraint is more flexible and favorable than the PIP constraint in the context 
of transmission over fading channels. Let Qav denotes the average interference power limit tolerated by 
PU-RX, then the AIP constraint can be written as, E[gop{go, gi)] < Qav 

The following optimal power allocation problem that minimizes the outage probability of SU in a 
narrowband spectrum sharing with one PU, under both a long term average transmit power (ATP) constraint 
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at SU-TX and an AIP constraint at the PU-RX, was considered in 

min Pr{i log(l + 5-1^(5-0, fi-i)) < ^o} 
p{go,gij>o ^ 

s.t. E[p{go,gi)] < Pav, E[gop{go, gi)] < Qav (1) 

where Pav is the maximum average transmit power at SU-TX. 

With the assumption that perfect CSI of both go and gi is available at the SU-TX, the optimal power 
allocation scheme for Problem ([B is given by flU: p*{go,gi) = ^ when \*j + ^*jgQ < ^, and otherwise, 
where c = e^''° — 1, and Aj, /Xj are the optimal nonnegative Lagrange multipliers associated with the ATP 
constraint and the AIP constraint, respectively, which can be obtained by solving X*j(E[p{go, gi)] — Pav) = 
and fi}{E[gop{go, gi)] - Qav) = . 

However, the assumption of full CSI at the SU-TX (especially that of go) is usually unrealistic and difficult 
to implement in practical systems, especially when this channel is not time-division duplex (TDD). In the 
next section, we are therefore interested in designing a power allocation strategy of the outage probability 
minimization Problem ([T]) based on quantized CSI at the SU-TX acquired via a no-delay and error-free 
feedback link with limited rate. 

III. Optimum Quantized power allocation (QPA) with imperfect gi and go at SU-TX 

A. Optimal QPA with limited rate feedback strategy 

As shown in Figdl following our earlier work [|T3l . we assume that there is a central controller termed as 
CR network manager who can obtain perfect information of go and gi, from PU-RX at the PU base station 
and SU-RX at the SU base station respectively, possibly over fibre-optic links, and then forward some 
appropriately quantized {go,gi) information to SU-TX through a finite-rate feedback link. For further 
details on the justification of resulting benefits due this assumption, see [13]. Under such a network 
modelling assumption, given B bits of feedback, a power codebook V= {pi, . . . , pi} of cardinality L = 2^, 
is designed offline purely on the basis of the statistics of go and gi information at the CR network 
manager. This codebook is also known a priori by both SU-TX and SU-RX for decoding purposes. Given 
a channel realization {go,gi), the CR network manager employs a deterministic mapping from the current 
instantaneous {go,gi) information to one of L integer indices (let Z{go,gi) denote the mapping, which 
partitions the vector space of {go, gi) into L regions TZi, . . . JZl, defined as I{go, gi) = j, if {go, gi) G 
TZj, j = 1, . . . , L), and then sends the corresponding index j = I{go, gi) to the SU-TX (and the SU-RX) 
via the feedback link. The SU-TX then uses the associated power codebook element (e.g., if the feedback 
signal is j, then pj will be used as the transmission power) to adapt its transmission strategy. 
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Remark 1: Note that the CR network manager could be assumed to be located at the SU-BS for the 
current setup and in this case, the PU-BS simply has to cooperate with the SU-BS by sending the real-time 
full CSI information of qq. However, for future generalization of our work to a multi-cell cognitive network 
scenario, we assume that the CR network manager is a separate entity, which can obtain information from 
multiple PU-BS and SU-BS if necessary. 

Define an indicator function Xj, j = 1, . . . , L, as Xj = 1 if | log(l + QiPj) < tq, and otherwise. Let 
Pr{TZj), E[»\Rj\ represent Pr{{go,gi) G TZj) and E[»\{go, gi) G TZj], respectively. Then the SU outage 
probability minimization problem with limited feedback can be formulated as 



L 

Pj>0, T^i Vj ^ 



L L 

S.t. J2E[pj\nj]Pr{n,)<Pay, J2^^9oPj\nj]Pr{nj) <Qay. (2) 

j=i j=i 

Thus the key problem to solve here is the joint optimization of the channel partition regions and the power 
codebook such that the outage probability of SU is minimized under the above constraints. 

The dual problem of ^ is expressed as, maxA>o, fi>o g{^,lA — ^Pav — j^Qav, where A,/i are the 
nonnegative Lagrange multipliers associated with the ATP and AIP constraints in Problem and the 
Lagrange dual function g{X, fx) is defined as 

L 

g{X,fi)= min Vi?[X, + (A + /i^7o)p,|7^,]Pr(7^,) (3) 

Pj>0, Tlj, Vj ^ 

The procedure we use to solve the above dual problem is: 
Step 1: With fixed values of A and ji, find the optimal solution (power codebook and quantization regions) 

for the Lagrange dual function ([3]). 
Step 2: Find the optimal A and /i by solving the dual problem using subgradient search method, i.e, updating 
A, fi until convergence using 



A'+i = [V - a\Pa. -J2E[Pj\n,]Pr{n,))]^, 

L 

^'+1 = _ /3'(g^, - ^ E[goPj\n,]Prin,))] + , (4) 

where / is the iteration number, a', are positive scalar step sizes for the l-th iteration satisfying 
'^i = < ^ similarly for and = max(a;, 0). 

Remark 2: A general method to solve Step 1 is to employ a simulation-based optimization algorithm 
called Simultaneous Perturbation Stochastic Approximation (SPSA) algorithm (for a step-by-step guide to 



implementation of SPSA, see IfTTl ). where one can use the objective function of Problem ([3]) as the loss 
function and the optimal power codebook elements for each channel partition are obtained via a randomized 
stochastic gradient search technique. Note that due to the presence of the indicator function and no explicit 
expression being available for the outage probability with quantized power allocation, we can't directly 
exploit the Generalized Lloyd Algorithm (GLA) with a Lagrangian distortion, as we used in [13J, to solve 
Problem ([3]). SPSA uses a simulation-based method to compute the loss function and then estimates the 
gradient from a number of loss function values computed by randomly perturbing the power codebook. 
Note that SPSA results in a local minimum (similar to GLA), but is computationally highly complex and 
the convergence time is also quite long. 

Due to the high computational complexity of SPSA and its long convergence time to solve Problem 
([3]), we will next derive a low-complexity approach for solving Problem (|3]). However, due to the original 
problem Q not being convex with respect to the power codebook elements, the optimal solution we can 
obtain is also locally optimal. 

Let V = {pi, . . . , pl}, where pi > ■ ■ ■ > Pl > 0, and the corresponding channel partitioning TZi, ... 
denote an optimal solution to Problem Let p{Z{go, gi)) represent the mapping from instantaneous 
{%^9i) information to the allocated power level. We can then obtain the following result: 

Lemma 1: Let {vi, . . . ,vl} denote the optimum quantization thresholds on the gi axis (0 < vi < ■ ■ ■ < 
vl) and {si, . . . , sl-i} indicate the optimum quantization thresholds on the g^ axis (0 < si < • • ■ < 
Then we have Vj, j = 1, . . . , L — 1, p{I{go, gi)) = pj, if vj < gi < fj+i, < go < sj and pi otherwise, 
where vj = = 1, . . . , L, and for Vj, j = 1, . . . , L — 1, when /i > 0, Sj = ^(^p Lp^-j — \^ while when 
/i = 0, = oo, then condition < (^o < boils down to A < ^ ^^^ ■ The region 71^ includes two parts 
: the set Uli = {{go, 9i) ■ vj < gi < Vj+i,go > Sj,yj = 0, . . . ,L - 1} with so = 0,Vo = and the set 
TIl2 = {(90,91) '■ 9i > vl,9o > 0}. The entire set TZli is in outage. 

Proof: See Appendix |Al 

Remark 3: When /i > 0, which implies that the AIP constraint is active, from Lemma [H the optimum 
partition regions possess a stepwise structure, as shown in FiglU When fi = 0, i.e, the AIP constraint is 
inactive and only ATP constraint is active (we must have A > 0), Problem Q becomes a scalar quantization 
problem involving quantizing gi only, and Lemma [U reduces to : p{X{gi)) = pj, if Vj < gi < fj+i, Vj, j = 
1, . . . , L — 1, and pi otherwise, where A < — , Vj = 1, 2, . . . , L — 1 and the two sub-regions of TZl 

Pj PL 

become T^li = {(71 : < (71 < fi} and 1^12 = {91 '■ 9i ^ vl}, and TZn is in outage. Note that in this 
case we must have Qav > Pav, due to Qav > Yl^=iE[goPj\Tlj]Pr{nj) = Yl^^i E[pj\nj]Pr{nj) = Pay, 
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where the last equality follows from the fact the E[gQ\TZj] = E[gQ] = 1 since TZj is formed purely based 
on the values of gi, which is independent of gQ. Note also that one can easily prove the converse, that 
when Qav > Pav, one must have /i = 0. 

From Lemma [H (due to the fading channels being independently exponentially distributed with unity 
mean) Problem (|2]) becomes, 

L-l 

min P^l, = 1 - e-''^ + ^(e"'^^ - e-''^+^)e-^^ 

L-l 

s.t. PL + Y^iPi - PL){e-^^ - e-'^^+Oll - e~'0 < Pav 

L-l 

VL + - VL){e~^' - e-^^+^)(l - e-^^(l + s,)) < Qa, (5) 

i=i 

where P^^^ denotes the outage probability with B = \0g2L bits feedback QPA, vj = = 1,...,L 
and for Vj, j = 1, . . . , L — 1, when > 0, sj = ^(^p ^^p^) — whereas when = 0, sj = 00. Although 
the above optimization problem may be verified to be non-convex, we can employ the KKT necessary 
conditions to find a local minimum for Problem ([5]). Taking the partial derivative of first order of the 
Lagrangian of Problem ([5]) over pj,j = 1, . . . , L — 1, and setting it to zero, we can obtain 

(e-. - e-''^+0[A(l - + /^(l - e~'^{l + s,))] = e'^^ ^[f{pj-i) - f{p,)], I < J < L - 1; (6) 

where f{po) = 1 and f{pj) = {pj — + — e~'^^)), l<j<L — 1.© also can be rewritten as 
j = l,...,L-l, 

Vj - ln(l - 

Equating the partial derivative of the Lagrangian function of Problem (|5]) over p^ to zero gives, 

L-l 

^(e-. _ e~^^^^)[Xil - + ' e-'^{l + s,))] + e''"^ ^f{pL-i) = A + /x. (8) 
i=i 

Optimal values of A and jj. can be determined by solving 

L-l 

i=i 

L-l 

I^IPl + Y^^Pj - PL){e-^^ - e--^+^){l - e-^^(l + s,)) - Q^,] = (9) 

Thus, for fixed values A and ji, we need to solve the L equations given by ([7]), ([8]) to obtain the power 
codebook. Given pi and pi, from ^ we can successively compute P2, • • • ,Pl-i, and then we can jointly 



solve the equation (|7]) with j = L — 1 and equation ([8]) numerically for pi and p^. The optimal value of A 
and /i can be obtained by solving ^ with a subgradient method, i,e. by updating A and /i until convergence 
using (HI). One can thus repeat the above two steps (i.e, given A and fi find the optimal power levels, and 
then using the resulting optimal power levels update A and fx) iteratively until a satisfactory convergence 
criterion is met. This procedure can be formally summarized as: 

a) First, if Pav < Qav, we must have = 0, A > 0. Starting with an arbitrary positive initial value for 
A, solve ©, ([8]) to obtain a power codebook {pi, . . . ,Pl}, and then use this codebook to update A by 
dH). Repeat these steps until convergence and the final codebook will be an optimal power codebook 
for Problem 

b) If Pav > Qav, we must have > by contradiction (since if ;U = 0, we must have Pav < Qav)- 
Let A = 0, then solving KKT conditions gives an optimal value of /i and corresponding power 
codebook {pi, . . . ,p£,}. With this codebook, if Xl^=i — ^av, then it is an optimal 
power codebook for Problem Q. Otherwise we must have A > too, in which case, starting with 
arbitrary positive initial values for A and /i, obtain the corresponding power codebook {pi, . . . ,pl}, 
and then update A and fj. by dH). Repeat these steps until convergence and the final codebook will 
be an optimal power codebook for Problem dS]). 

B. Suboptimal QPA Algorithm 

When the number of feedback bits B (or alternatively, L) goes to infinity, we can obtain the following 
Lemma that allows us to obtain a suboptimal but computationally efficient quantized power allocation 
algorithm for large but finite L. 

Lemma 2: ImiL^^pL = 
Proof: See appendix |Bl 

Remark 4: Lemma [21 shows that regardless of whether /i > or = 0, with high rate quantization, the 
power level for the last region TZc approaches zero, which also implies the following as L — )■ oo: 

1) The non-outage part of TZc, given by TZc2^ disappears gradually. In other words, TZc — Tlci- Thus, 
when L — )• oo, TZc becomes the outage region with zero power assigned to it. 

2) When fi > 0, the quantization thresholds on the axis sj — )■ s'j (where s'j = ^ — which gives 
Vj = cA + cfis'j, and it means all the points given by coordinates {s'j, Vj) lie on the line of gi = cX + cfigo. 
Therefore, as L — )• oo, the stepwise shape of the structure in yU > case (i.e, the boundary between 
non-outage and outage regions) approaches the straight line gi = cX + cngo, which is consistent with the 
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full CSI-based power allocation result in dU. 

Thus, when L is large, applying Lemma [2l (i.e, pi — )■ 0) to Problem ©, the above L KKT conditions 
Q and dH]) can be simplified into L — 1 equations: 

(e-. - e-'^^+i)[A(l - e-''^) + - e'^'^il + s'^))] 

= e-^'^^lPj-iiX + /^(l - e-^^-O) - Pj{X + /^(l - e-^^))], Vj, j = 1, . . . , L - 1 (10) 

where when fi > 0, the quantization thresholds on the go axis are given by s'j = ~ ^» ■^g = 0, and 
Po = x+^is' ' while when /i = 0, s'j = oo, Sq = 0, and po = ^. (flOl ) can be also written as 

c 

Pi+i = — , j = !,•••, ^-2; 

^b,_i(A+M(l-e ^.-i))-p^(A+Ml-e 

- ln(l - ^ ^ ^ ) 



A(l-e J)+A»(l-e 

A(i-e-^-i) + /i(l-e-'^-i(l + .U)) 
^[pi_2(A + - e-^^-)) - p^_i(A + /i(l - e-^i-i))] 



1 (11) 



Thus, for given values of A and fi, starting with a specific value of pi, we can successively compute 
P2, ■ ■ ■ ,Pl-i using the first equation of (fTT]) (recall that vj = j-). Then the second equation in (fTTI) 
becomes an equation in pi only, which can be solved easily using a suitable nonlinear equation solver. We 
call this suboptimal QPA algorithm as 'Zero Power in Outage Region Approximation' (ZPiORA), which 
is applicable to the case of a large number of feedback bits, where the exact definition of "large" will be 
dependent on the system parameters. Through simulation studies, we will illustrate that for our choice of 
system parameters, ZPiORA performs well even for as low as 5 = 2 bits of feedback. 
Alternative suboptimal algorithms: For comparison purposes, we also propose two alternative suboptimal 
algorithms described below: 
(1) The first suboptimal algorithm is based on an equal average power per (quantized) region (EPPR) 
approximation algorithm, proposed in [.18,1 in a non-cognitive or typical primary network setting 
for an outage minimization problem with only an ATP constraint. More specifically, by applying 
the mean value theorem (similar to IfTSlI ) into the KKT conditions © with j = 2, . . . , L — 1, 
we can easily obtain that Pjie^""^ — e"''-'+^)[A(l — e""*^) + — e~'*j(l + Sj))] ^ pj-ii^e'""^-^ — 
e""'^)[A(l - e-^i-^) + - e"*J-i(l + sj^i))], j = 2, . . . , L - 2. Adding the two equations of © 
together and applying ([8]), we have J2j=i Pji^'''^ ~ e"'"^+^)[A(l — e^''^) + — e~'*^(l + Sj))] = 
APa.+/iQa.-e-"^^(pL_i-pL)(A+^(l-e-^^-0). Since e-''^^(pL-i -PL)(A+/x(l-e-^^-0) can 
be approximated as PL-ii^^^'^'^ — e~^^)[A(l — e^"*^-!) + — 6""*^-^ (1 + Si_i))] by using the mean 
value theorem, we can obtain the following L (approximate) equations, namely pj (e~^J — e~^^+^)[A(l — 
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e-'^)+fi{l - e-'^{l + Sj))] ^ hPav+HQa^^ j = 1, . . . ,L - 1 and - Et=i (e^''^ - e-^^+0[A(l - 
e~^^) + — e^*^ (1 + Sj))]) _ Then one can jointly solve the above L equations and two 

other equations with j = 1 and ([8])) for A, /x, pj, Vj = 1, . . . , L. We call this suboptimal algorithm 
as the "Modified EPPR (MEPPR)" approximation algorithm. Obviously, ZPiORA is computationally 
much simpler than this method, especially when > 0. Furthermore, from simulations, when Pav 
or Qav is small, the performance of ZPiORA is always better than MEPPR. It is seen however that 
when both Pav and Qav are large, for a small number of feedback bits, MEPPR may outperform 
ZPiORA, whereas with a sufficiently large number of feedback bits, ZPiORA is a more accurate 
approximation due to Lemma [2] (when L is large, approaches zero, whereas MEPPR has pl > 
VL). See Section V for more details. Note that, an equal probability per region (excluding the outage 
region) approximation algorithm employed in [jTOll for scalar quantization can not be applied to our 
case (vector quantization), since it will increase the computational complexity even further. 
(2) The second algorithm is based on GLA with a sigmoid function approximation (GLASFA) method 
proposed by lfT9l . where the sigmoid function is used to approximate the indicator function in 
the Lagrange dual function ([3]). More specifically, given a random initial power codebook, we use 
the nearest neighbor condition of Lloyd's algorithm with a Lagrangian distortion d{{go, gi), j) = 
Xj + (A + figo)pj to generate the optimal partition regions [20J given by, TZj = {{go, gi) : Xj + (A + 
I^QojPj < Xi + {\ + fj.go)pi,yi 7^ j}, i,j = 1, . . . ,L. We then use the resulting optimal partition 
regions to update the power codebook by pj ^ argminp^>gi?[(T(/i;(| log(l + giPj) — r^)) + (A + 
^go)pj\Rj]Pr{Rj) for j = 1, . . . , L, where we use the approximation Xj ^ (t(A;(| log(l + giPj) — 
To)), o-(x) = being the sigmoid function where the coefficient k controls the sharpness of the 
approximation (for detailed guidelines on choosing k see [19|). The above two steps of GLA are 
repeated until convergence. Numerical results illustrate that ZPiORA significantly outperforms this 
suboptimal method. See Section V for more details. 

IV. Asymptotic outage behaviour of QPA under high resolution quantization 

In this section, we derive a number of asymptotic expressions for the SU outage probability when the 
number of feedback bits approaches infinity. To this end, we first derive some useful properties regarding 
the quantizer structure at high rate quantization: 

Lemma 3: As the number of quantization regions L — oo, we can obtain the following result: with 
fi > 0, the optimum quantization thresholds on the go axis satisfy 
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s[ -s'q^ s'2-s[^ ■■■ ^ - s'^_2, where s'j = ^ - ^, j = 1, L - 1 and s'q = 0. With /i = 0, the 
optimum quantization thresholds on the gi axis satisfy ^ ^ ^ ■ ■ ■ ^ tif^' where Vj = j-, j = 1, . . . , L — 1 
and here t>o = cA. 

Proof: See Appendix O 
Lemma 4: In the high rate quantization regime, as L — ^ oo, we have 



3=1 i=i 



(12) 



where when /i > 0, = ^ — ^, whereas when /i = 0, = oo, and (fT2l) simplifies to ce ~ '^^=1 
with f,- = — . 

Proof: See Appendix iDl 

With Lemma [3] and Lemma IH the main result of this section can be obtained in the following Theorem. 



Theorem 1: The asymptotic SU outage probability for a large number of feedback bits is given as, 

-i(i+^) 

PLt ~ 1 — e'" / [I — (1 — e~^) ^ 1 ] (for ;U > 0) where a is a constant satisfying 

CM J 

1-e E 

CUI -a(l+^) c{u*fY _a(l+^) 1 

^[(A} + /.})(l--^(l-e 1^^^., (1-e -/\l + a(l + —)))]. (13) 



We also have limi^oo^ott = 1 - e"''^/[l - ^'\^^ ^ ]• For /i = 0, P„tt ~ 1 - e"'^/(^+^\ where /3 is a 
tant given by e^""^/ ?^ A} 
Proof: See Appendix |El 



constant given by e "^^^ ?^ AjPai)^^-^-'-- In this case we also have limi_>oo Pj^t = 1 — e "^'^ 



V. Numerical Results 

In this section, we will examine the outage probability performance of the SU in a narrowband spectrum 
sharing system with the proposed power allocation strategies via numerical simulations. All the channels 
involved are assumed to be independent and undergo identical Rayleigh fading, i.e, channel power gain 
and gx are independent and identically exponentially distributed with unity mean. The required transmission 
rate is taken to be ro = 0.25 nats per channel use. 

Fig. [3] displays the SU outage probability performance of the suboptimal algorithm ZPiORA versus Pav 
with feedback bits B = {1,2}, under Qav = —5 dB and Qav = dB respectively, and compares these 
results with the corresponding outage performance of the suboptimal method MEPPR and the optimal QPA. 
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As observed from Fig. [3l when Qav = —5 dB, with B fixed, the outage performances of ZPiORA and 
corresponding optimal QPA almost overlap with each other. When Qav = dB and Pav < —5 dB, with the 
same number of feedback bits, the outage performances of these two methods are still indistinguishable; 
and with Pav > —5 dB, the outage performance gap between ZPiORA and corresponding optimal QPA is 
decreasing with increasing B. For example, with 1 bit feedback, at Pav = 10 dB, the outage gap between 
ZPiORA and optimal QPA is 0.0347, but with 2 bits of feedback, the outage performance of these two 
methods are very close to each other, which agrees with Lemma |2] that ZPiORA is a near-optimal algorithm 
for large number of feedback bits. Now we look at the performance comparison between ZPiORA and 
MEPPR. As illustrated in Fig. [3l when Pav or Qav is small, with B bits feedback, the performance of 
ZPiORA is always better than MEPPR. This is attributed to the fact that when Pav or Qav is small, it can 
be easily verified that is close to zero, but MEPPR always uses > 0. However, when both Pav and 
Qav are large (e.g. Pav > dB and Qav = dB), for 1 bit feedback case, MEPPR outperforms ZPiORA 
and performs very close to the optimal QPA, whereas with a sufficiently large number of feedback bits (in 
fact, with more than just 2 bits of feedback), ZPiORA is a more accurate approximation due to Lemma [2l 
These results confirm the ZPiORA is a better option for a large number of feedback bits, not to mention 
that ZPiORA is much simpler to implement than MEPPR. 

In addition. Fig. |4] compares the outage performance of ZPiORA with another suboptimal method 
(GLASFA) with Qav = — 5 dB. We can easily observe that with a fixed number of feedback bits (2 
bits or 4 bits), ZPiORA always outperforms GLASFA. And ZPiORA is also substantially faster than 
GLASFA. For example, with fixed A and n and 4 bits of feedback (Qav = — 5dB, Pav = 10 dB), when 
implemented in MATLAB (version 7.11.0.584 (R2010b)) on a AMD Quad-Core processor (CPU P940 
with a clock speed of 1.70 GHz and a memory of 4 GB), it was seen that GLASFA (with 100,000 training 
samples, starting k = 20 and increasing it by a factor of 1.5 at each step which finally converged at about 
k = 768.8672) took approximately 299.442522 seconds (different initial guesses of the power codebook 
may result in different convergence time). In comparison, ZPiORA took only 0.006237 seconds to achieve 
comparable levels of accuracy. These results further confirm the efficiency of ZPiORA. 

Fig. |5] illustrates the outage performance of SU with optimal QPA strategy versus Pav with feedback 
bits B = {2,4,6}, under Qav = —5 dB and Qav = dB respectively, and studies the effect of increasing 
the number of feedback bits on the outage performance. For comparison, we also plot the corresponding 
SU outage performance with full CSI case. Since ZPiORA is an efficient suboptimal method for large 
number of feedback bits, we employ ZPiORA to obtain the outage performance instead of using optimal 
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QPA for B = 6 bits. First, it can be easily observed that all the outage curves decrease gradually as Pav 
increases until Pav reaches a certain threshold, when the outage probability attains a floor. This is due 
to the fact that in the high Pav regime, the AIP constraint dominates. For a fixed number of feedback 
bits, the higher Qav is, the smaller the resultant outage probability is, since higher Qav means PU can 
tolerate more interference. Fig. [5] also illustrates that for fixed Qav, introducing one extra bit of feedback 
substantially reduces the outage gap between QPA and the perfect CSI case. To be specific, for Qav = 
dB and Pav = 10 dB, with 2 bits, 4 bits and 6 bits of feedback, the outage gaps with the full CSI case are 
approximately 0.1083, 0.0249 and 0.006979 respectively. And for any Qav, only 6 bits of feedback seem 
to result in an SU outage performance very close to that with full CSI case. 

Figure [6] compares the asymptotic outage performance derived in Theorem \T\ and the optimal QPA 
performance B = {4, 6, 8} under Qav = dB. It is clearly observed that increasing number of feedback 
bits substantially shrinks the outage performance gap between the asymptotic outage approximation and 
the corresponding optimal QPA performance. For instance, with 4, 6, 8 bits of feedback at Pav = 10 
dB, the outage gap between the asymptotic outage approximation and the corresponding optimal QPA are 
around 0.0325, 0.00618, 0.000168 respectively. These results confirm that the derived asymptotic outage 
expressions in Theorem [T] are highly accurate for B > 8 bits of feedback. In addition. Figure |7] depicts 
the asymptotic outage probability behavior of SU versus the number of quantization level L at Qav = 
dB, Pav = 10 dB, and compares the result with the corresponding full CSI performance. It can be seen 
from Figure |7] that the outage decreases as the number of quantization level L increases, however, as L 
increases beyond a certain number (L > 2^, i.e, B > 8 bits), the outage probability curve starts to saturate 
and approaches the full CSI performance. This further confirms that only a small number of feedback bits 
is enough to obtain an outage performance close to the perfect CSI-based performance. 

VI. Conclusions and extensions 

In this paper, we designed optimal power allocation algorithms for secondary outage probability min- 
imization with quantized CSI information for a narrowband spectrum sharing cognitive radio framework 
under an ATP constraint at SU-TX and an AIP constraint at PU-RX. We prove that the optimal channel 
partition structure has a "stepwise" pattern based on which an efficient optimal power codebook design 
algorithm is provided. In the case of a large number of feedback bits, we derive a novel low-complexity 
suboptimal algorithm ZPiORA which is seen to outperform alternative suboptimal algorithms based on 
approximations used in the existing literature. We also derive explicit expressions for asymptotic behavior 
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of the SU outage probability for a large number of feedback bits. Although the presented optimal power 
codebook design methods result in locally optimal solutions (due to the non-convexity of the quantized 
power allocation problem), numerical results illustrate that only 6 bits of feedback result in SU outage 
performance very close to that obtained with full CSI at the SU transmitter. Future work will involve 
extending the results to more complex wideband spectrum sharing scenario along with consideration of 
other types of interference constraints at the PU receiver. 

Appendix 

A. Proof of Lemma U} 

We use an analysis method similar to [1211 to prove our problem's optimal quantizer structure. Let V = 
{pi, . . . ,pl}, where pi > ■ ■ ■ > Pl > 0, and the corresponding channel partitioning 71 = {TZi, . . . , TZl} 
denote the optimal solution to the optimization problem ((21), and p{go,gi) = Pj, if ido^gi) ^ T^j- 

Let n* = {(go^gi) : vj < gi < t;,+i,0 < g^ < s,}, j = 1, . . . , L - 1 and 7^2 = 7^2l U Tll^ = 
{{go^gi) : Vj < gi< vj+i, go > Sj^Mj = 0, 1, . . . , L - 1} U {{go^gi) : gi > VL,go > 0}, where so = and 
Vq = 0. We assume that the set TZ* \ TZj is a non-empty set, where \ is the set subtraction operation (i.e, 
if (fi'O; gi) thsii (do, gi) ^T^] but ((70, gi) ^ TZj). Then, the set TZ* \ TZj can be partitioned into 

two subsets SJ = (7^* \ TZj) n (U^7^fc) and 5+ = (7^* \ TZj) n (uj^^j^^TZk). In what follows, we denote 
the empty set by 0. 

(1): We will show that Sr = 0, Vj = 1, . . . , L. 

(a) : When j = 1, it is obvious that = 0. When 1 < j < L, if Sj' ^ 0, then we can always reassign 
the set Sf into region TZj without changing the overall outage probability. This is due to the fact that 
within the set Sj' G TZ*, we have vj < gi < fj+i resulting in ^ log(l + giPj) > tq, and the power level in 
(u{~^7?.fc) satisfies pk > pj. Thus Sj' is never in outage. However, the new assignment can achieve a lower 
Lagrange dual function (LDF) in ([3]), due to g'{X,fi) — g{X,fJ,) = E[{X + fJ,go){pj — Pk)\Sf]Pr{Sj') < 0, 
where g'{\,n) denotes the LDF with the new assignment, which contradicts the optimality of the solution 

r,n. 

(b) When j = L, if 5^ ^ 0, we can again reassign the set into region TZl- 1) If some part of is in the 
set {{go, gi) ■ < gi < vi, go > 0} of 7^2i' we have | log(l+5'ipi) < ro, which implies that this part of 

is always in outage. Therefore, this reassignment for this part of S£ will not change the outage probability 
but will decrease the LDF due to the power level pl in TZi is the lowest. 2) For any j (j = 1, . . . , L — 1), 
if some part of (denoted by "S^^") exists in the set {{go,gi) : Vj < gi < t'j+i, go > Sj} of TZ*li, we 
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have l\og{l + giPj) > tq, ^ log(l + 5(1^^+1) < tq and {\ + figQ){Pj -Pl) > 1- And given 5"^^^ C {U^lJ^TZk), 
let the power level for 5*^^ be pk (where k could be any value from — 1}). Reassigning this 

part of set Sj^ into region IZl will reduce the value of the LDF, since if k < j (implying pk > Pj), 
^'(A, yu) - ^(A, /i) = E[l + pl{\ + /i^o) - Pk{\ + ^igi))\Slp]Pr{Slp) < and if A; > j (implying pk < pj), 
g'{X,IJ,)~g{X,fi) = E[1+pl{X + figo) - 1 -pk{X + lJ-go)\Slp]Pr{S;[;p) < 0. 3) If some part of belongs 
to the set 7^22' similar to (a), we can show that the new partition for this part of does not change the 
overall outage probability and meanwhile reduces the value of the LDF. These all contradict optimality. 
(2): We will now show that the set 5^ = 0, j = 1, . . . , L. When j = L, it's straightforward that S'l = 0. 
When j < L, we assume that Sj' 7^ 0. Within the set Sj' E TZ*, we have Vj < gi < Vj+i, implying 
|log(l + giPj+i) < To, or in other words, Sj' G (U^^j^-^^TZk) is in outage. We can reallocate the set Sj' 
into region TZj. This reassignment not only lowers the outage probability (Sj' with pj will not be in outage) 
but also lowers the value of the LDF, given by g'{X, /i) —g{X, fi) = E[{X + figo){pj ~Pk) — l\S^]Pr{S^) < 
E[{X + fJ,go){pj —Pl) — l\Sj']Pr{Sj') < 0, due to go < sj = ^^^pLp^) — \- This also contradicts optimality. 

Therefore, we have 7^* \ 7^J■ = 0, Vj = 1, . . . , L, i.e, 7^* C Hj, Vj = 1, . . . , L. Since = 
the whole space of {g^, gi) = \J^^{IZj, and TZ* C 7^^, Vj, we can obtain that TZ* = TZj, Vj = 1, . . . , L. 

B. Proof of Lemma |2} 

We assume that lim^^^ooPL 7^ 0. Let 5 = limL-^oo Pl > 0. From the KKT condition ([8]), we have 

L-l 

e-^-^{pL-i -Pl){X + - e--0) = (A + f^){Po'ut + e"^^) + /^^^^"'^ " e-''^^')^-'^ 

>(A + ^)(P,^„, + e-^-) (14) 

Let P/„( denote the outage probability with full CSI at SU-TX, then we have P^^t > PLt and lim^-^oo PLt = 
P/„^. Taking the limit L — > 00 on both sides of (fT4l) . we have 

lim e-^'^^ipL-i - Pl){X + /i(l - e-^^-)) > (A + fi){PL + e't) ^ (15) 

L^oo PJ^ 

Given pi >■■■> pl > 0, it is clear that the sequence {pj}, j = 1, 2, . . . , L is a monotonically decreasing 
sequence bounded below, therefore it must converge to its greatest-lower bound 5, as L — )■ 00. Therefore, it 
can be easily shown that for an arbitrarily small e > 0, we can always find a sufficiently large L such that 
Pl-1-Pl < e. Thus, as L ^ 00, (pl-i-Pl) 0, which implies when /j, > 0, sl^i = ^ oo- 

This implies that 

Mm e-''^^{pL^i-pL){X + fi{l-e-'^-^)) = e-l^iX + fi) lim (p^-i - p,.) = 0. (16) 

L-)-oo p"^ 0^ L^oo 
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which is in contradiction with (fT5l ). Thus, we must have MmL^ao Pl = 0. 



C. Proof of Lemma \3\ 

As L — 7- oo, from Lemma[2l we have pl — )• 0. Applying it to Problem ([5]), we have the KKT conditions 
as dH). 

1) /i > 0: From = ^ - ^, we have pj = and we also have po = j^^^. Applying it to (HO]), 

the right hand side (RHS) of equation (flOl) becomes, 



RHS = ,-.4|A + Ml-e-^-) _ AH-Ml-e- 



^ ^ ^ A+/j(l-e _ A+/^(l-e J ) 



From the mean value theorem (MVT), we have 

A+^t(l-e~"i-i) _ A+^(l-e~°J) 

where s' E [s'j_i, s'j). As the number of feedback bits B = logg L — i- oo, the length of quantization interval 
on gQ axis [s'j_i, s'j),j = 1, . . . , L — 1 approaches zero IfTSl . Hence (fTSl ) becomes. 



A+At(l-e "J-I) _ A+At(l-e ) 

A+/iS^_j A+/is^ — /i 



— [A(l-e-^)+Ml-e-^Hl + <))] (19) 



Applying ^ to ([H]), we have RHS ^ e-''^cfi{s'j - s;_i)[A(l - e""^) + - e-"^(l + s^)]. Similarly, 
as L — 7- oo, we also have the length of quantization interval on the gi axis [vj,Vj+i),j = 1, ... ,L — 2 
approaches zero, thus from MVT, e^'"^ — e^^^+^ ~ e'""^ (^i+i Thus the left hand side (LHS) of equation 
(fTOl) can be approximated as, LHS ^ e~'"^{vj^i — Vj)[\{l — e^^'i) + — e"*^(l + s'j)). Hence, we have 
Vj = 1, . . . , L - 2, Vj+i - Vj ^ c/i(s^. - from which we get s'j^^ - s'j ^ s'j - Vj = 1, . . . , L - 2, 

namely, s'i_i — s'^_2 ^ ■ ■ ■ ^ s[ — s'q, since vj = cA + c/is^. 

2) /i = 0: In this case, we have s'j = oo,j = 1, . . . , L — 1. Thus (flOl) becomes e"*'^ — e~*'j+i = 
e~''-' -%(j9j_i — pj), where j = 1, . . . , L — 1 and po = j, which can be rewritten as — (e~^J — e~^^+'^) = 
;jj:^e~^j(fj — fj-i), where = ^ = cA. Applying MVT into as before, we have ^e~'"^ (f j+i — Vj) ~ 
^e-^^ (i;j-t;j_i),Vj = 1, L-2, which yields ^ ^ ^.Vj = 1, L-2, namely, ^ ~ ■ ■ ■ ~ ^• 

This completes the proof for Lemma [3l 
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D. Proof of Lemma 

As L — i- oo, from Lemma [2l we have p/, — 0. Adding the two equations of ^ together and applying 
Pl — 0, we have 

L-l 

J^p,(e-^^ - e-^^+0[A(l - e-^^) + - e-^^l + 4))] = APa. + y^Qa. (20) 
The KKT conditions (flOl) can be rewritten as 

p,(e-''^ - e-'^^+OfAll - e-^'^) + - e^^Hl + ^;))] = - ^^_,)IZ^fcil^ZM (21) 

where f'{pj) = Pj(A+yu(l — e"*^)). As mentioned before, when L — t- oo, we have the length of quantization 
interval on the gi axis [vj-i, Vj),j = 2, . . . , L — 1 approaching zero. Hence we also have the length of the 
interval [pj-i,Pj),j = 2, . . . , L — 1 approaching zero, since Vj = ^. Thus from MVT, we have 

fip,.,) - f'{p,) ^ _ ^ _ ^ ^22) 

Pj-i - Pj 

Applying (|22l) into (|2TI) . we can obtain, Vj = 2, . . . , L — 1 

p,{e-'^ - e-^^+0[A(l - e-''^) + fi{l - e-^^(l + s'^))] ^ 

p,_i(e-^-i - e-^0[A(l - e"'-0 + /^(l - e-^^-^(l + (23) 
Then applying the result of (l23l) into (l20l) . we can have j = 1, . . . , L — 1 



which gives, 

g(e-. _ e-^-)[A(l - e-^^) + /i(l - e'^Hl + ^;))] - ^-^- + ^/^- -• 
i=i i=i 
This completes the proof for Lemma HI 



(25) 



E. Proof of Theorem [7]' 

1) ;U > 0: From Lemma [H we can easily obtain, s'j ^ js[, j- = X + fis'j ^ X + j^s\, and Vj = j- 
cX + jc^is'^,^] = 1, . . . , L - 1, Let z = E^Ti^le""^ - e-^^+i)[A(l - e~''^) + - e-"^(l + s'^))], which 
implies that < z < A + u. Then from Lemma HI we have -r^ T^f i^ — ~ z' , where z' = ,„ f ^ 
and < z' < , „ . Using the above results, we get s\ ^ Hl^ziA) = | where ci = Let 
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a = cfid = 2{z' — X)c, then s[ ^ Since < z' < xp^+^q~' ha\e lim2,_j>oo f = 0. From the 
definition of z above, we have 

L-l L-l 

^ e-'\{X + fi)e-T - (1 - e-*)(A + " (1 - e-^)M je-^^^M)] 

i=i i=i 

L-l L-l 

-cA r 

\ A -i- H. \H ^ II H ^ll>\-t-//. I > H " 1^ II i^, 

where 6 = a + Ls'^ = a{l + ^) and we also have lim^^oo x = 0- Since Yj^=i ^^^^ = 1 — 1' ^^d 

1— e 

A_M,_ A 



- e -^[(A + ^)e-i - (1 - e~T)(X + /i) e"^':^ - (1 - e't)^ J] je"^-^] (26) 



EL-l . _jA e'"L"''^(Lei -e''-L+l) j^zi u 
je = ^ 5 '-^ ij2M becomes 

(l-e~"r)2 



b . 



z ^ e-^^[(A + ;x)(l - (1 - e-T)- -) - (1 - g-r) — ^ -I L] (27) 



Since Um^^oo f = and limi^oo |; = 0, we have 1 — e ^ ~ f and 1 — e i ^ |-. And when L — t- oo, 
we approach the full CSI scenario, thus implying A ^ A^,/i ^ /i-'^. Using these results in (|T7l) . we have 

-cA* 



a , , , b 



z ^ e ^'}[{X} + f,}){l - -(1 - e-^)) - ^((1 - - e-') - be-')] 



^ e~^'}[{X} + /i*)(l - ^(1 - e-^)) - ^(1 - e-\l + 6))] (28) 

Since z = {XPav + fJ'Qav)z' = {XPav + fiQav){X + ^ {X}Pav + fJ'jQ av) {X} + fj, wc Can obtain a from 
the following approximation: 

iX}Pa,+fi}Qa.)iX} + ^y'*f 

CU*f -a(l+^) c{u*fY -a(l+^) 1 

From (|29l ). with given Pa^, and Qat;, a is a constant. Then when L is large, 

L-l L-l 

Pout ~ 1 - e""^ + 5Z(e~"^- - e-^^+^e^'^- ^ 1 - e-'^[e-i - (1 - 6-^)5^6^^'^] 

-a(l+-t^) 

= 1 - - (1 - e-*)-^^] ^ 1 - e-'^^ni " (1 - e'^) .,,,^J (30) 



l-e 



-a(l+^) 

and limi->oo ^ott = 1 - [1 " ] 



2) /i = 0: Let y = ^ = ^ > 1, then again, from Lemma [3l we can get for j = 1, . . . , L — 1, 
cA?/-'. From Lemma H we have e""^-^^ ^ ^T^H^jZlv^ = ^i^i?- With x = y - 1, we have, 

x'(L-l) 
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Vj ^ cA?/-'. i^rom Lemma |4' ^""''^ — M^v^L^ 

g-cA(i+x) ^ APa„(l + a:) ^^^!^-^'^ - Now, suppose limi^ooa;^ = oo. Since (1 + x)^-^ > I + [L - l)x 



|(L — 2)(L — we have liiiii^^oo ^"^^L~i)a~^ ^ limL-i>oo 1 + |(-^ — 2)x = oo. Then taking the limit as 
L — )■ oo, we have limi_>oo e~'^'^^^+^^ = cxd, which contradicts linii^oo e"'^'*'^^"^^'^ < 1, thus we must have 
liuiL^ooxL = < P < oo (where /3 is a constant), implying as L -> oo, x — )■ |. Applying this result, 
we get c^'^^'^^'^l) ^ XPavi^ + f )^~^77 — r~"- After taking the limit as L — )■ oo on both sides of above 
equation, we have e~'^^f ~ X*j:Pav^-f^, from which one can solve for /3 approximately. Note that in the 
above approximation, we have used lim2,^oo(l + f )^~^ = and when L is large, A ~ Aj. Therefore, 
when L is large, Pj;^ = 1 - e^^^ = 1 - 6-^=^(1+^) ^ 1 - e"'^/(^+i\ limL_,oo PLt = 1 " e~^^^'. 
This completes the proof for Theorem \T\ 
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Fig. 1. System model for narrowband spectrum sharing scenario with limited rate feedback 
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Fig. 2. The 'stepwise structure' of optimum quantization regions for fi > case 




Fig. 3. Outage probability performance comparison between ZPiORA, MEPPR and optimal QPA 
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Fig. 4. Outage probability performance comparison between ZPiORA and other possible suboptimal algorithm : GLASFA 
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Fig. 5. Effect of increasing feedback bits on outage performance of SU 




Fig. 6. Comparison between asymptotic outage performance and QPA performance with Qav = OdB 
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. Q - Asymptotic outage behaviour at O^^-O dB, P^^^l dB 
fuli CSi outage at Q -0 dB, P -1 dB 
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Fig. 7. Asymptotic outage beliaviour versus the number of quantization level L 
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